On the relation between the boson-fermion model and RVB state 



Evgucny Kochetov 

Bogolyubov Theoretical laboratory, Joint Institute for Nuclear Research, 141980 Dubna, Russia 



(N 
O 
O 
(N 

Oh' 
< 

m 



I 



o 

> 
O 

O 
(N 
O 



O 

o 



X 



Marcin Mierzejewski 
Institute of Physics, University of Silesia, 40-007 Katowice, Poland 



We show that the resonating valence bond (RVB) state 
of the t-J Hamiltonian can be described within an effective 
boson-fermion model. We derive a boson-fermion Ifamilto- 
nian which describes the lowest order coupling between holons 
and spinon waves. The RVB order parameter determines a 
hybridization between the both types of excitations. The 
effective Hamiltonian is investigated within the lowest or- 
der self-consistent conserving diagrammatic approximation. 
We compare our approach with the phenomenological boson- 
fermion model (BFM) that is believed to describe the pseu- 
dogap phase. 
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The origin of the pseudogap is one of the most im- 
portant current problems in high temperature super- 
conductors. This normal state gap is seen by various 
techniques, including angle-resolved photoemission spec- 
troscopy (|-||], NMR infrared j^, and transport 
measurements. Although we now possess consider- 
able experimental data, the understanding of this phe- 
nomenon is still far from complete. There are a number 
of microscopic scenarios, but none of them is commonly 
accepted. The opening of the pseudogap, and its evolu- 
tion into the true superconducting gap, may, however, be 
described by a phenomenological boson-fermion model 
|p|-|T2t. This BFM is usually understood as a system of 
localized tightly bound electron pairs (bosons) which hy- 
bridize with itinerant electrons. This model customarily 
neglects the hard core effects of the bound electron pairs 
|p|-pl| and assumes the standard commutation relations 
for the boson operators. We may then identify three 
characteristic temperatures pO[ | which correspond to: (i) 
the opening of the pseudogap in the fermion excitation 
spectrum, (ii) the appearance of itinerant bosons, (Hi) 
the condensation of these bosons and the concurrent on- 
set of superconductivity. Due to the phenomenological 
character of the BFM, the microscopic origin of the pair 
correlations is not specified. Moreover, it is not straight- 
forward to select an appropriate set of model parameters 
— in particular the total number of fermions and bosons. 
In order to answer these questions we need to find a re- 
lationship between the BFM and a microscopic approach 
p3[ . The aim of this work is to provide such a link by 
demonstrating a close similarity between the BFM and 
the t-J model with a mean-field RVB state. 

As it is known, the RVB state is characterized 
by three important features: (i) spin-charge separation. 



(ii) short range antiferromagnetic (AF) correlations, and 
{Hi) holon condensation. Owing to («), the RVB pair- 
ing amplitude may develop at a much higher tempera- 
ture than the superconducting transition temperature, at 
which the phase coherence between the electron pairs is 
established p5[ | . In the RVB scenario there is a phase 
with the preformed, uncoherent, electron pairs with an 
amplitude Arvb that serves as the RVB order parameter. 
In this phase the RVB electron pairs are to be thought 
of as spinon singlets which are coupled to holon excita- 
tions, this resulting in a boson-fermion (BF) interaction. 
It is therefore natural to suggest that it is the RVB order 
parameter that controls the pseudogap. This is our basic 
assumption. 

First we demonstrate that the RVB state can be de- 
scribed by a holon-spinon interaction. Next, applica- 
tion of the linear spinon-wave approximation leads to a 
boson-fermion model where bosons are doubly charged 
with respect to spinless fermions. 

To begin, we express the t-J model in terms of the 
Hubbard operators ATf^ = cj^(l — ni__o-). 



H.C. + J 



where Cia- annihilates an electron at site i with spin 

projection a = ±, rn^ = cJo-Cio-, = Sergio- and 
denotes a summation over the nearest neighbor 

nonrepeated bonds. The electron-spin operator, Qj = 



2 L-iaa' 

and r 



X^Ta^iX^l" , enters the magnetic part of 



t-J 



(T^,r^,r'^) stand for the Pauli matrices. The 
operators Q fulfil the SU(2) algebra commutation rela- 
tions, but do not form a true representation, since the 
second Casimir operator is not a c-number. 

Note that the X'^" appear as fermionic-like operators, 
while the X'^'^ = X'^^X^'^ correspond to bosonic de- 
grees of freedom. The crucial t-J model requirement 
of no double occupancy (^la^ia ^ 1 follows automati- 
cally: by definition X°° + J2a = I- Unlike the slave- 
particle representations, this restriction imposes an ad- 
ditional constraint on dynamical degrees of freedom and 
seems to pose a severe problem in a mean-field treat- 
ment. Note further that the t-J Hamiltonian possesses 
two global U(l) symmetries: \Jn^{1) and Uq3(1). These 
correspond to the conservation of the total electron num- 
ber — '^ii^i'^ — TT'ii a-nd to the total spin 
projection Q3 = ^ J2ii'^t^ ~ -^i )' respectively. 
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A conventional Hartree-Fock decoupling is applied to 
the magnetic part of the t- J Hamiltonian 



Hj := J 



QiQj - ^niUj 



to yield 



ihj) 



|A '2 



where &J^. = ^ - , is the valence 

bond " singlet" pair creation operator and Ay is the RVB 
order parameter defined on each bond between the near- 
est neighbor sites. A.^ / breaks the global UjVe(l) 
symmetry, though it does not directly result in supercon- 
ductivity. It instead indicates the onset of the electron 
spin-singlet formation. There are many possible mean 
field states in RVB theory, among which the extended 
s-phase is selected here because of its simplicity. 

To proceed, we employ the Feynman path-integral 
representation for the partition function. We use the 
geometric quantization approach based on the Poisson 
realization of the SU(2|1) superalgebra on a classical 
phase space — the supersphere CP^'^ In this ap- 

proach there is an isomorphic map X — > X'^^ of the 
SU(2|1) generators into classical functions on the su- 
persphere which preserves the Lie superalgebra struc- 
ture. This map can be realized explicitly by exhibit- 
ing covariant symbols (known as the momentum maps) 
of the Hubbard operators in terms of the SU(2|1) co- 
herent states, @: = {z,^\X\z,0- Then, 
= -i(l ~ \z\^)w, {Q+y^ = zw, (g-)^' = zw, 

{x"-y^ = -z^w, ix°+Y^ = -^w, {x+'^Y^ = -Cw, 
{x-^y' = -z^w, {x++Y^ + {x—y' =_ (i + zz)w, 

{X°°Y^ = £_£_w, where w := (1 + \z\^ + ^^)-\ Here 
the even z and odd ^ Grassmann variables parameter- 
ize superspace SU(2|1)/U(1|1) = CP^^Ii), which is the 
A'^ = 1 supersymmetric extension of the Riemannian 
sphere CP^. These functions form a SU(2|1) superal- 
gebra under the Poisson superbrackets on CP^'^. For an 
arbitrary fermionic Hamiltonian with the requirement of 
no double occupancy, the relevant Berezin path-integral 
representation of the partition function is equivalent to 
the SU(2|1) phase-space path integral 

Z = j -D^su(2|i)exp (^J{z,S.\-^ - H\z,S_)dt^ , 



'DAisu(2|i) 



dzdzd^d^ 
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-CO' 



where the Hamiltonian is expressed in terms of the X- 
operators. Despite a rather complicated appearence, 
the SU(2|1) path-integral representation for the t-J 
model is convenient since (i) it automatically incorpo- 
rates the requirement of no double occupancy, {X^^^^ + 



J^ai'^'^'^T'' = independetly of any approximations be- 
ing applied to z and ^ fields; (ii) it furnishes a convenient 
set of independent dynamical variables {z,£,) £ CP^'^, 
which can directly be related to the charge and spin de- 
grees of freedom. 

The coordinates z and ^ can be related to the SU(2) 
spinon and U(l) holon degrees of freedom in the follow- 
ing manner: We make two successive changes of vari- 
ables z z\J\ -I- ^ ^-^/l + |zp to decouple the 
SU(2|1) measure, -Dmsu(2|i) ^Msu(2)^Mu(i), where 
D^su(2) = dzdz I 'Ki(\^\z^Y stands for the SU(2) invari- 
ant measure on CP^, while -D/iu(i) = dS^dS, denotes the 
Berezin integration over odd Grassmann variables [ pTf . 

The t-J partition function becomes 



Zt-j - 



(^^^ n ^'^su(2) / n exp[A-j] 



(1) 



with corresponding to the U(l) charged spinless 
fermion degrees of freedom (holons), and Zi representing 
pure SU(2) spins (spinous). 

The classical Hamiltonian now reads 

^ArvB 



- /i' ^(1 - U^) - [253(z„ z.) -t- ■ (2) 

i i 

We have dropped the constant term. We have also ex- 
plicitly introduced a chemical potential term, — /i'A^g, as 
well as the Langrange multiplier A which controls the 
magnitude of the total electron magnetic moment. 

The SU(2) algebra is realized on CP^ in terms of 

the momentum maps Sf = — ^ Yqqfp, (S*^)^' = 



^{Zj,Z^) 



{Zj Zi) 



appears as a spinon-singlet amplitude, $y = 

where ^'s3=±i/2(-2^) is a spinon wave function in the 
basis of the SU(2) coherent states, \z) = (1 + 
|z|2)-i/2exp{z5_}|.S3 - 1/2). Namely, ^^.^^^./^{z) = 
(S3 = ±1/21^), 5(3) |S3) = S3IS3) and the SU{2) 
coherent-state inner product is 



{z,\zj) 



(1 -f ZiZj) 
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so that = 1 - \{zi\zj)\'^. 

In general, the Hamiltonian (H) cannot be thought of as 
a covariant symbol of a Fermi-Bose type operator, except 
for a particular case of the linear spinon-wave (LSW) ap- 
proximation. This approximation amounts to expanding 
the action (Q), as well as a measure factor in the path 
integral, in powers of |zp up to the leading order. In the 
paramagnetic phase (Qff^) = 0, and the spinon fluctu- 
ations become bounded by the condition (l^il^) < 1/2. 
This, to some extent, justifies the LSW approximation. 
Within this approximation the partition function, takes 
on the form 

Zbf = / d^jd^j exp[^BF], 

j 

-4bf = 2 ^ X V^^^ ^ ^^^^ + ^^^^ ~ ^^^^ ^ ^ 
where the classical Hamiltonian reads: 
^b'f = E ^^^^ ~ E ^^^^ (^"^ - ^"^) + H-c. 

- ^i' E(i - - ^E + 2^^-^^ + ■ 

i i 

It straightforward to see that the above partion function 
represents the following boson-fermion system: 



' 1 df 

BP I "''J 



Hbf ^ -tY,M] -'^^Y.f^^M 

(id) (id) 

+ ABj2bih-^,j2i2bib, + fln), 

i i 



H.c. 



(3) 



The parameters of the boson-fermion system are ex- 
pressed in terms of the parameters of the original t~J 
Hamiltonian and the RVB state: v — JArvb/4; = 
2fi'; fi — ^' + X. The Uq3(1) symmetry of eq. (||) un- 
der global gauge transformations Zi — > e^^^Zi, e^^^i 
leads to conservation of the total charge J2ii'^blbi + fj fi). 
Note, that this global U(l) symmetry characterizes the 
standard BFM and indicates that bosons are doubly 
charged with respect to fermions. The BFM discussed 
in Refs. [p|-pl| is refered to as the standard approach. 
One should also keep in mind that this symmetry is an 
intimate feature of the RVB state, and occurs indepen- 
dently of the LSW aproximation. 

The standard BFM was investigated with the help of 
various methods, e.g.: the lowest order self-consistent 
conserving diagrammatic approximation 10|, the dy- 
namical mean-field approximation and the contin- 
uous unitary transformation [ pT| . Unlike the standard 
BFM, the Hamiltonian obtained in the present study in- 
volves spinless fermions. Therefore decay and recombina- 
tion between fermions and bosons involves two neighbor- 
ing lattice sites. In the present case, bosons with momen- 
tum q = do not participate in these processes, in the 



standard BFM these bosons are of crucial importance. 
Because of these differences, we cannot directly adopt 
the numerical results obtained for the standard BFM. 
Therefore, we have investigated a one dimensional (ID) 
case within the lowest order self-consistent conserving 
diagrammatic approximation and compared our results 
to those presented in Refs. ||9|,p^. 

In the the lowest order self-consistent conserving dia- 
grammatic approximation one obtains a system of equa- 
tions for the fermion, GfO^,oj), and boson, G'B(k, w). 
Green's function. We refer to Ref. |]lO| for the details. 
This set of equations has been solved iteratively slightly 
above the real axis, i.e. for oj — Keuj + iO.Olt. We have 
used grids of 501 momenta in the interval [0, tt] , and 2001 
frequencies in the interval [— 5t, 5i]. We have also fixed 
the bosonic level Ab = —2.41. 
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FIG. 1. The fermionic density of states calculated at the 
temperature kT — 0.02i for different values of the bo- 
son-fermion coupling, u) = corresponds the Fermi level. 
The inset shows the average bosonic occupation number as a 
function of momentum. 

Fig. 1 shows the fermionic density of states, p(tf), cal- 
culated for different values of the RVB order parame- 
ter. The boson-fermion coupling enhances the density 
of states just above the Fermi level and reduces p for 
much larger energies. Effectively, the local minimum in 
the density of states occurs at the Fermi level. Modifica- 
tion of the density of states takes place also for a weak 
boson-fermion coupling, when the local minimum at the 
Fermi level is almost not visible. Therefore, despite the 
absence of the gap the system cannot be considered as a 
composition of noninteracting fermions and bosons. The 
most important difference between the standard BFM 
and the present study shows up in the average boson oc- 
cupation number ub (see the insets in Figs. 1 and 2). 
Here, the boson-fermion coupling enhances the number 
of bosons with momentum close to tt, whereas in the stan- 
dard BFM the bosons with momentum q = play the 
dominating role. In the present case the momentum de- 
pendence of rt_B can be attributed to the AF short range 
correlations. However, bosons are uniformly distributed 
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in the real space and there is no antiferromagnetic order 
(in agreement with the assumption of the RVB state). 
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FIG. 2. The fermionic density of states calculated for 
(i) non interacting particles, (ii) standard BFM (fermions 
with spin) and (iii) obtained in the present study (spinless 
fermions). The inset shows the average boson occupation 
number as a function of momentum in the cases (ii) and (iii). 

In Fig. 2 we compare p{uj) and ub obtained within 
the standard BFM and the present study. In the case of 
spinless fermions riB(k) is relatively flat, when compared 
to the standard approach. It suggests that bosons do not 
condensate at the state with some particular momentum. 
However, investigation of this problem is not possible 
within a one-dimensional numerical study. One should 
also keep in mind, that in our approach (which concerns 
only the mean field RVB state) condensation of bosons 
does not represent the superconducting phase transition. 
Similarly to the standard BFM, the local minimum in p, 
that occurs at the Fermi level, deepens with the decreas- 
ing temperature and disappears when the temperature is 
sufficiently high (see Fig. 3). 

The difference between the present approach and the 
standard BFM shows also up in the fermionic spectral 
functions (see Fig. 4). Here, the most important broad- 
ening of the spectral functions takes place above and be- 
low the Fermi surface. In the standard case ^ the spec- 
tral functions have a three-peaked structures, which are 
mostly visible at the Fermi surface. It may be attributed 
to the free-particle-like motion of the bosons with small 
wave vectors. However, in the present approach these 
bosons are not involved in the decay and recombination 
processes. 
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FIG. 3. The same as in Fig. 1 but calculated at different 
temperatures for v = O.lSt. 

To conclude, starting from an RVB state for the t-J 
Hamiltonian, we have derived a connection between the 
t-J model and the BFM. We have expressed the param- 
eters of the effective Hamiltonian in terms of those of the 
t~J model and the RVB state. The resulting model dif- 
fers from the standard one in that the fermionic degrees 
of freedom correspond to spinless quasiparticles. How- 
ever, as in the standard BFM, a local minimum in the 
fermionic density of states occurs close to the Fermi level. 
The minimum deepens with the decreasing of tempera- 
ture and vanishes when the temperature is sufficiently 
high. The present study is restricted only to the mean- 
field RVB state, and does not attempt to account for the 
superconducting phase transition. It is, however, free 
from the ambiguities which occur in the phenomenologi- 
cal BFM when tries to select an appropriate set of model 
parameters. 
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FIG. 4. Spectral function for fermions calculated for 
wave-vectors, which are in the vicinity of the Fermi surface. 
We have used v = 0.15t and kT = 0.02t. 
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